base_dir = "/Users/pol/owncloud/HiWi_folder/02_running_projects/02_MA_thesis/01_Data/Preprocessed/"
pt_path = paste0(base_dir, "PitchTiers/08_manual_corrected/")
Fujisaki_path = paste0(base_dir, "Fujisaki/")
library(contouR)
Compute slopes by drawing a linear regression line. Let’s first get an intuition by looking at the first 5 slopes.
files = list.files(pt_path)
files = files[grepl(".PitchTier", files)]
compute_regression_slope(files[1:5], pt_path, plot_only=TRUE)
## Loading required namespace: hqmisc
Now compute the slopes for all files:
features = compute_regression_slope(files, pt_path)
head(features)
## name a b reg_RMSE naive_intercept naive_slope
## 1 DF_ANG_VX1 26.14186 -9.7273180 2.572510 23.25355 -9.93577139
## 2 DF_ANG_VX10c 17.91798 0.1759461 3.313159 17.19068 2.13530953
## 3 DF_ANG_VX11a 22.49966 -3.4898652 3.705590 18.66072 -2.47602240
## 4 DF_ANG_VX12 21.70759 -2.1948611 2.689266 15.86251 0.05407044
## 5 DF_ANG_VX13a 28.59098 -10.0375849 3.931166 24.21989 -11.33068258
## 6 DF_ANG_VX14 20.33091 -4.5572098 4.454911 20.07695 -8.69039835
## naive_RMSE min_max_slope min_max_RMSE
## 1 3.948228 -17.84009 4.188648
## 2 3.485802 -57.95575 28.774414
## 3 4.762231 -25.15022 14.607542
## 4 4.776461 -32.97042 18.349630
## 5 6.796655 -17.42301 5.193256
## 6 5.479758 -17.08930 8.409380
Let’s check the correlations. The regression intercept correlates strongly with the naive intercept. There is also a smaller positive correlation between naive and regression slope. However, it is not correlated with the slope between global minimum and maximum.
library(corrplot)
## corrplot 0.84 loaded
interesting_columns = 2:9
interesting_columns = names(features)[interesting_columns]
cor_df = features[, interesting_columns]
print(corrplot(cor(na.omit(cor_df)), method="number", type = "upper"))
## a b reg_RMSE naive_intercept
## a 1.000000000 -0.41068879 0.46425991 0.86460664
## b -0.410688789 1.00000000 -0.03463125 -0.18848927
## reg_RMSE 0.464259908 -0.03463125 1.00000000 0.42803043
## naive_intercept 0.864606643 -0.18848927 0.42803043 1.00000000
## naive_slope -0.487538388 0.53115991 -0.41912562 -0.54656440
## naive_RMSE 0.418484349 0.01567463 0.86820015 0.33379923
## min_max_slope 0.003059545 0.09329694 -0.06747011 -0.05259831
## min_max_RMSE 0.110679211 0.09607382 0.12344647 0.05688564
## naive_slope naive_RMSE min_max_slope min_max_RMSE
## a -0.48753839 0.41848435 0.003059545 0.11067921
## b 0.53115991 0.01567463 0.093296944 0.09607382
## reg_RMSE -0.41912562 0.86820015 -0.067470114 0.12344647
## naive_intercept -0.54656440 0.33379923 -0.052598310 0.05688564
## naive_slope 1.00000000 -0.40522122 0.188722278 0.07816047
## naive_RMSE -0.40522122 1.00000000 -0.031021720 0.13827723
## min_max_slope 0.18872228 -0.03102172 1.000000000 0.92176783
## min_max_RMSE 0.07816047 0.13827723 0.921767829 1.00000000
features$filename = features$name
features = add_meta_data(features)
interesting_columns = c("a", "naive_slope", "b", "naive_intercept")
for (col in interesting_columns){
print(significance_test(na.omit(features), col))
}
## Loading required package: ggplot2
## Warning: package 'ggplot2' was built under R version 3.5.2
## Loading required package: magrittr
From these optimal hyper parameter combinations we can compute features on the phrase command
top_scores = read.csv("top_scores.csv")
phrase_features = compute_phrase_features(top_scores, Fujisaki_path)
names(phrase_features) = c("name", "num_phrases", "Ap", "phrase_T0", "phrase_reg_slope", "phrase_reg_RMSE", "phrase_naive_slope", "phrase_naive_RMSE")
head(phrase_features)
## name num_phrases Ap phrase_T0 phrase_reg_slope
## 1 DF_ANG_VX1 1 0.73900 -0.46 -0.003836904
## 2 DF_ANG_VX10c 2 0.07395 NA NA
## 3 DF_ANG_VX11a 1 0.25830 -0.47 -0.001377021
## 4 DF_ANG_VX12 1 0.28640 -0.46 -0.001097925
## 5 DF_ANG_VX13a 1 1.18280 -0.34 -0.005092311
## 6 DF_ANG_VX14 1 0.65570 -0.47 -0.003469269
## phrase_reg_RMSE phrase_naive_slope phrase_naive_RMSE
## 1 0.009014474 -0.0034409529 0.022443360
## 2 NA NA NA
## 3 0.034902102 -0.0015639914 0.079482718
## 4 0.005779488 -0.0009455769 0.009563485
## 5 0.038778261 -0.0040847909 0.067799577
## 6 0.022504489 -0.0031715768 0.038879622
features = combine_features(features, phrase_features, key = 'name')
We see that both slope computations (and also the RMSE) are highly correlated. Also we find a negative correlation between the amplitude and the slope
corrplot(cor(na.omit(phrase_features[, 3:8])), method="circle", type = "upper")
An visualize the features in boxplots:
phrase_features$filename = phrase_features$name
interesting_columns = names(phrase_features)[c(2:5, 7)]
for (col in interesting_columns){
print(significance_test(na.omit(phrase_features), col))
}
## Warning: Computation failed in `stat_compare_means()`:
## data are essentially constant
We can also devide the contour in it’s ICCs. Let’s first give some examples for the ICC separation:
compute_EAC_ICCs(files[1:5], pt_path, plot_only=TRUE)
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
We can now extract the features from the ICC as Mirjana Rajkovic, Jovicic, Grozdic, Zdravkovic and Subotic (2018) propose
feature = compute_EAC_ICCs(files, pt_path)
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(0). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(0). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(0). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in compute_EAC_ICCs(files, pt_path): No ICC detected
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in compute_EAC_ICCs(files, pt_path): No ICC detected
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(0). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(0). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(0). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(0). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(0). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(3). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(1). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
## Warning in warn0("ICC size to small(", end_idx - start_idx, "). Current ICC skipped"): ICC size to small(2). Current ICC skipped
head(feature)
## name num_ICCs num_pos_ICCs num_neg_ICCs mean_delta_F0
## 1 DF_ANG_VX1 3 0 3 6.763778
## 2 DF_ANG_VX10c 4 1 3 5.867933
## 3 DF_ANG_VX11a 6 1 5 6.336774
## 4 DF_ANG_VX12 6 1 5 6.107005
## 5 DF_ANG_VX13a 3 0 3 9.458832
## 6 DF_ANG_VX14 4 0 4 8.457135
## sd_delta_F0 mean_tc sd_tc mean_delta_tc sd_delta_tc EAC
## 1 1.120784 0.2214286 0.14146645 0.2008929 0.12259159 NA
## 2 4.720117 0.2544643 0.20634854 0.2057143 0.20932654 0.003093019
## 3 4.255035 0.1767857 0.14047569 0.1545918 0.14104075 0.070351529
## 4 1.638377 0.1732143 0.08474903 0.1836735 0.08216504 0.053752638
## 5 7.077855 0.4535714 0.42077734 0.3535714 0.39753706 NA
## 6 6.395935 0.2303571 0.26455103 0.1885714 0.24742810 NA
features = combine_features(features, feature, key = 'name')
We can have a look at the correlations:
# Many
feature$filename = feature$name
cor_df = feature[, names(feature)[2:10]]
corrplot(cor(na.omit(cor_df)), method="circle", type = "upper")
cor_df = feature[, names(feature)[2:11]]
corrplot(cor(na.omit(cor_df)), method="circle", type = "upper")
An visualize the features in boxplots:
interesting_columns = names(feature)[2:11]
for (col in interesting_columns){
print(significance_test(na.omit(feature), col))
}
Last but not least, we have to save the features in a CSV for later use:
write.csv(features, "slope_features.csv", row.names = FALSE)